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Carbon nanotubes (CNTs) have large intrinsic carrier mobility due to weak acoustic phonon scat- 
tering. However, unlike two-dimensional metal-oxide-semiconductor field effect transistors (MOS- 
FETs), substrate surface polar phonon (SPP) scattering has a dramatic effect on the CNTFET 
mobility, due to the reduced vertical dimensions of the latter. We find that for the Van der Waals 
distance between CNT and an Si02 substrate, the low-field mobility at room temperature is reduced 
by almost an order of magnitude depending on the tube diameter. We predict additional experi- 
mental signatures of the SPP mechanism in dependence of the mobility on density, temperature, 
tube diameter, and CNT - substrate separation. 



Semiconducting carbon nanotubes (CNTs) show 
promise for technological applications in electronics and 
optoelectronics primarily due to weak acoustic phonon 
carrier scattering and a direct bandgap [ij . The low-field 
mobilities in CNTs have been measured by several groups 
0, S 0] • At low biases only acoustic phonons contribute 
to the scattering in CNTs and the scattering is weak 
0, 0, 0]. The CNT diameter d scaling of the low-field 
mobility fj, was predicted Q to be ^ (x cP and confirmed 
experimentally [3|. However, the magnitude of the exper- 
imental mobility was found to be an order of magnitude 
smaller than the theoretical prediction. Recently it was 
shown that surface polar phonon (SPP) scattering can be 
strong in CNTs lying on polar substrates such as Si02 
[1] . While the SPP scattering is of lesser importance in 
the conventional two-dimensional MOSFETs 0, as we 
show below it is much more prominent in CNTs due to 
the much smaller vertical dimension of the devices given 
by the CNT diameter. 




FIG. 1: Schematics of the SPP scattering in CNTs. 



In this Letter we evaluate the transport properties 
of CNTs due to the extrinsic scattering mechanism by 
surface phonons in the insulating polar substrates as 
schematically shown in Fig. [TJ The surface phonons on 
polar substrates produce a surface electric field coupled 
to the electrons on the nearby CNTs. We find that SPP 
scattering lowers dramatically the low-field mobility, even 
at room temperature. What is more surprising, the SSP 
scattering can increase the drift velocity under high bias 
conditions, especially at high doping levels and elevated 
CNT phonon temperatures. 

The interaction potential is given by [lol |: 
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where R and r are substrate phonon and CNT electron 
coordinates, respectively. The strength of the SPP scat- 
tering is given by the polarization field P{R) at the sur- 
face associated with a particular excited phonon mode; 
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where at^ is a creation operator of a surface phonon 
mode v with two-dimensional wavevector q. The vector 
normal to the solid (at x < 0) is x and the unit vector 
along q is q. The normalization surface area A drops out 
from the final result. The magnitude of the polarization 
field is controlled by the Frohlich coupling F,^: 
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where hjjso.v is a surface phonon energy and eo and Coo 
are the low- and high-frequency dielectric constants of the 
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polar substrate ll|. In the following we evaluate scat- 
tering in CNTs on Si02 surfaces, which have five phonon 
modes huJso,^=i...5 = 50,62,100,149,149 meV, with 
Frohlich couplings ^2=1.. .5 = 0.042, 0.38, 0.069, 1.08, 1.08 
mcV respectively [12[. 

The CNT electron-polar surface phonon matrix ele- 
ment has been derived in the envelope wavefuncton ap- 
proximation in 
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where 4'^ and '^k+q are electron wavefunctions in CNT 
with momenta k and k + q charaterized by the quantized 
angular momentum m and a continuous z component 
along the tube axis, d is a tube diameter, and h is the 
height from the CNT surface to the insulating substrate, 
see Fig. [T] The length of the tube L determines the k- 
point sampling along the tube axis Aq^ = 27r/L. The 
functions / and K are the Besscl functions of the first 
and second order, respectively. It was shown that SPP 
scattering is strongest when the azimuthal angular mo- 
mentum is conserved. Am = 00]. Therefore, to analyze 
the results of the calculations it is instructive to con- 
sider the asymptotic form of the scattering potential for 
Am = and large qz'. 



47re2F2 



L 



qzd^/2TT{d + 2h) 



(5) 



This form of the SPP coupling also works well for the 
intermediate values of q^. 

The elcctron-CNT phonon scattering is described using 
the Su-Schriefer-Heeger model as in ref. Q with the cou- 
pling constant g = 5.3 eV/A. The electron band struc- 
ture is described by the 7r-orbital tight-binding model 
with to = 3 eV. 

We study the effect of the SPP scattering on transport 
in the diffusive regime. We calculate the mobility by 
solving the steady-state multiband Boltzmann transport 
equation (BTE) in the presence of an electric field. The 
temperature dependent scattering rates Wk.k+q between 
electronic states k and k+q entering in the BTE equation 
are evaluated according to: 

+ (1 + (5(efc - Efc+q - /ia;_q^.)] (6) 

where the phonon occupation number is given by 
the Bose-Einstein distribution. The phonon index v runs 
over all phonon modes involved in scattering, i.e. six 
CNT phonon modes and five surface phonon modes in 
Si02, the latter are treated here as being dispersiveless. 

The inverse mobility at low-field and low density is 
shown in Fig. [2] as a function of temperature. At low 
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FIG. 2: Temperature dependence of the low-field, low density 
inverse mobility in (19,0) tube [d = 1.5 nm) on a log-log scale. 
The red circles are in the presence of both CNT phonons and 
the SPP scattering, the cyan circles with CNT phonons and 
one surface phonon at hw^ = 62 meV, blue circles with the 
CNT phonon scattering alone. The black solid curve is a 
fit to Eq. (0 with acNT = 0.0013 and aspp = 13.4. The 
low temperature mobility is influenced by the CNT acoustic 
phonon scattering and gives linear temperature dependence 
shown by the green dashed line. The high temperature inverse 
mobility can be approximated by the power-low fits shown by 
the green dashed (T^'®) and dashed-dotted lines (T^'^^). 



temperatures, the inverse mobility is linear in tempera- 
ture, because only CNT acoustic phonon scattering dom- 
inates transport in that regime. At temperatures above 
about T=100 K, SPP scattering is activated and starts 
to dominate transport. The temperature dependence of 
the low-field mobility becomes super linear, due to the 
activated nature of the SPP scattering and can well be 
approximated by the following formula: 
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where acNT is a contribution from the acoustic CNT 
phonons and aspp is proportional to the SPP scattering 
rate, ns is the occupation number of SPP with ener- 
gies hLtj2 = 62 meV, huji^s = 149 meV, and P ~ (F| + 
F^)/F2 — 5.6 is the ratio of the electron-polar phonon 
coupling strengths fixed by the ratio of the Frohlich cou- 
plings. Since the occupation number of the highest en- 
ergy SPP phonon of 149 meV is much smaller than the 
huj2 =62 meV mode, in a zeroth order approximation 
the mobility below room temperature can be calculated 
through a coupling only to the 62 meV SPP phonon 
mode. We note however that this approximation is not 
exact and at least the two lowest energy modes hiui and 
hu!2 have to be included to get an adequate mobility at 
room temperatures. The single mode 11102 calculations 
shown in Fig. [5] (cyan color) overestimate mobility by 
50% at room temperature. 

Despite the fact that the two low energy SPP phonons 
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at room temperature have small occupations of ub = 
0.15 and 0.09 respectively, they do dominate the mobility, 
because of much stronger coupling than that to CNT 
acoustic phonons. The mobility at room temperature 
degrades by an order of magnitude in the presence of 
the SPP mechanism, depending on the tube diameter 
and its distance from the substrate surface h. The latter 
depends ultimately on the roughness of the surface. The 
low-field mobility dependence on h can be rationalized 
by Matthiessen's rule, with the SPP scattering rate given 
by Eq. (O and CNT phonon scattering described by the 
mobility ficNT in isolated (suspended) CNT: 



fi = 
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where iJLsi02 is mobility on polar surface substrates with 
the minimum distance determined by the Van der Waals 
interaction = 0.35 nm. The electron momentum 
transfer along the tube axis in the scattering of a surface 
phonon is given by hvpg^ ~ \Jfujj2{2/^ + fujj2), where 
2A « 0.9/(i eV nm is a bandgap and vp ~ 10^ cm/s 
is the Fermi velocity. We find that the ratio of the low 
field mobilities at room temperature can be as large as a 
factor of 8.6 in 2.5 nm diameter tubes. 




as found in ref. [1, Q and Fig. 3 in Supporting Informa- 
tion. At higher temperatures, SPP scattering contributes 
to mobility with a weaker diameter dependence than that 
for CNT phonons. This results in smaller a in Eq. © at 
elevated temperatures. Therefore, the relative contribu- 
tion to the total mobility of the SPP phonons in Fig. O 
is larger in larger diameter tubes, i.e Hcnt / IJ-Si02 ^ dP'^ 
at T=300 K and ^icNT/^J-Si02 oc d° '^ at T=500 K. 



FIG. 3: Dependence of the low-field, low-density mobility on 
the inverse distance from the polar substrate at room tem- 
perature on a logarithmic scale. The mobility is reduced by a 
factor of f^CNx/ fJ'Si02 = 5.3, 6.5, 7.5, and 8.6 correspondingly 
in (13, 0) - red, (19,0) - blue, (25, 0) - green, and (31, 0) - cyan 
tubes. The solid curves are calculations according to Eq. (jSj 
with (7z=0.38, 0.32, 0.28, and 0.25 nm"\ 

The diameter dependence of the low-field, low-density 
mobility shows a power law behavior: 



H{d) = Ad", 



(9) 



with power a being temperature dependent as shown in 
Fig.dl At low temperatures, only acoustic phonons in the 
CNT contribute to the scattering and the power is a « 2 




1.4 1.6 1.8 2.0 2.2 2.4 
d (nm) 



FIG. 4: Low-field, low-density mobility as a function of tube 
diameter at three temperatures on a log-log plot for all chi- 
rality tubes. The black solid lines are best fits to Eq. ^ with 
{A,a) =(2.6,2.18), (0.23,1.67), and (0.069,1.47) for T = 100 
K - blue, T = 300 - green, and T = 500 - red respectively. 

Although it is natural to analyze low-density mobility 
theoretically, it is difficult to measure it experimentally 
because the device is close to the off-state with a high 
resistance. In the case of only CNT phonon scattering, 
as the device is electrostatically or chemically doped with 
carriers, the mobility was shown to increase due to the 
smaller density of states available for scattering as the 
Fermi level moves away from the bottom of the conduc- 
tion or valance bands [13[ . As the Fermi level approaches 
the second band, a new channel for scattering opens up 
and the mobility drops showing a non-monotonic behav- 
ior. In the presence of the SPP scattering, the density de- 
pendence of the mobility shows an additional minimum, 
when the Fermi level approaches the lowest energy SPP 
phonon energy ^012 with the strongest coupling, as shown 
in Fig. [5^. The size of the mobility modulation depends 
on the temperature T, which determines the Pauli block- 
ing factors (1 — gk) ~ e^^'^"""/'^^-^, where gk is the 
distribution function from the BTE solution and Epermi 
is the quasi-Fermi level. At low temperatures T < 100 K, 
the Pauli blocking principle does not allow any SPP scat- 
tering even if the energy conservation law can be satisfied 
E Fermi > hjJi. At higher temperatures the mobility min- 
imum broadens out with increasing temperature. As the 
density further increases, the mobility increases until it 
reaches the bottom of the second band, showing a second 
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minimum associated with the higher subband scattering 
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FIG. 5: (a) Carrier density dependence of the low-field mo- 
bility in the same tube normalized to its value at zero den- 
sity. The vertical lines shows energies of the two polar optical 
phonons and the second subband, at which scattering time is 
being reduced. The upper axis shows the energy of the Fermi 
levels for a given density at zero temperature, (b) Resistance 
per unit length in (19,0) tube {d — 1.5 nm) at different tem- 
peratures 30 K (black), 50 K (blue), 100 K (cyan), 150 K 
(green), 200 K (magenta), 300 K (red) on a log scale. 

We have calculated the dependence of the resistivity 
on carrier density for the same tube and the results are 
shown in Fig. [SJd. The resistivity (resistance per CNT 
tube length) is non-monotonic with a minimum corre- 
sponding to occupation of the second subband at low 
temperatures and a second minimum emerging at inter- 
mediate temperatures T sa 100 K due to the SPP scat- 
tering. The prediction of the appearance of a second 
minimum in the resistivity at higher temperatures could 
provide an experimental signature of the SPP mecha- 
nism. On the other hand, a deviation of the temperature 
dependence of the low- field mobility from the 1/T law 
is a more subtle evidence for the SPP mechanism, since 
the position of the threshold voltage can be tempera- 
ture dependent due to the hysteresis [3, 15 1. In addi- 
tion, the maximum mobility associated with the onset of 
the scattering in the second subband shows a superlinear 
temperature dependence (/x oc T^'^) even for only CNT 
phonon scattering see Fig. 1 in Supporting Information. 
At the same time in metallic tubes, the density of states 
around the neutrality point is constant and, therefore, 



a much weaker density dependence of the mobility (or 
conductance) is expected, which makes the temperature 
dependence of the mobility to be a more suitable test for 
the SPP mechanism. Indeed, in metallic tubes, a non- 
linear terngerature dependence of the resistivity has been 
reported 
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FIG. 6: High field velocity (at F=5 V//im) as a function of 
carrier density in three tubes (13, 0) d=1.0 nm - red, (19, 0) 
d=1.5 nm - green, (25, 0) d=2.0 nm - blue. Solid curves are 
obtained with both CNT and Si02 phonon scattering, dashed 
curves with only CNT phonon scattering, and dotted curves 
with both CNT and Si02 phonon scattering, but carries re- 
stricted to the first subband. Polar phonons are kept at T = 
300 K, while CNT acoustic and optical phonons are allowed 
to heat up to (a) Tac = 300 K Top = 1000 K, (b) T= 1000 K 



and (c) Tac 



1000 K respectively. 



Finally, we discuss the transport properties of semi- 
conducting CNTs under the high bias conditions. At 
low-densities with only CNT phonon scattering, the drift 
velocity Vhigh was found to saturate with the field F at 
about half the Fermi velocity Vgat ~ 5x 10^ cm/s (for 
the choice of to=3 eV): ■y^/g^ = (^o^^)"^ + f^at H- As 
the density increases, the high field velocity drops and 
density dependence of the velocity does not show a sim- 
ple saturated form [l^]. Here we report high field ve- 
locity calculated at a constant F = 5 V/^m, which can 
be achieved in the experiments before the long channel 
CNT burns under high bias condition [H, [13, 111. The 
results are shown in Fig. [6l The reference dashed lines 
are calculations with only CNT phonon scattering. As 
the SPP scattering process is turned on, the high bias 
velocity drops at small densities, especially in smaller di- 
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ameter tubes, while it increases at larger densities and 
larger diameter tubes as shown in Fig. [6^. This is due 
to both the hyperboHc dispersion of the CNT bandstruc- 
ture [l^ and the multiband nature of transport at high 
biases. The SPP scattering is much stronger than the 
CNT optical phonon scattering and, therefore, the tail 
of distribution function or effective electronic temper- 
ature [i^l (Te eFXph) is expected to be smaller in 
the presence of the SPP mechanism, due to the shorter 
mean free path \ph- We find that at low densities, first 
subband dominates transport by using restricted calcula- 
tions with electrons occupying only the first subband as 
shown in Fig. [6] (dotted curves). As the density increases, 
the band velocity increases and so docs the drift velocity. 
In the multiband calculations, the higher subbands with 
higher effective masses and lower band velocities influ- 
ence transport, which lowers the drift velocity in Fig. [6] 
(solid curves) . In the presence of only CNT phonon scat- 
tering, even at low densities the high bias transport is 
influenced by higher subbands, which makes the quan- 
titative interpretation to be more difficult although it 
shows similar trends. 

It has been suggested that under high bias condi- 
tions the optical phonons can be out of equilibrium with 
the acoustic CNT phonons 17,[2l|. Recently, the non- 
equilibrium phonon population in CNTs under high bias 
conditions was directly observed experimentally [22] . To 
study the influence of the hot phonon effect on the trans- 
port properties, we calculate the drift velocity as a func- 
tion of density for three different tubes in Fig. ^p. In 
the case of only CNT scattering, we find a significant 
drop in the drift velocity with the elevation of the op- 
tical phonon temperature. This drop gives rise to the 
negative differential conductance as observed in the sus- 
pended CNTs [l7|. On the other hand, the magnitude 



of the drift velocity change in the presence of the SPP 
is much smaller, because SPP, which dominates trans- 
port, are kept at ambient temperature. The weaker ef- 
fect of the CNT phonons on transport properties leads 
to near disappearance of the negative differential con- 
ductance even when CNT optical phonon temperature 
is elevated Finally, we find that the drift velocity 

under high bias conditions in the presence of the SPP 
scattering is not effected by the acoustic CNT phonon 
temperature as is seen from Fig. |6};, which is essentially 
identical to Fig. [6b. 

In conclusion, we calculated the effect of the surface 
polar phonon scattering on transport in semiconducting 
CNTs. We find that at room temperature SPP scattering 
reduces the mobility by an order of magnitude from the 
intrinsic value determined by the CNT acoustic phonon 
scattering. The diameter dependence of the mobility in 
the presence of SPP can still be well approximated by a 
power low, but with a reduced temperature dependent 
value of the power as the relative role of the SPP scat- 
tering on transport increases with increasing tube diam- 



eter. At temperatures above about 100 K, we predict an 
additional minimum in the mobility-density dependence 
associated with the low energy SPP phonon scattering. 
This can serve as an experimental signature of the SPP 
mechanism. Under high bias conditions we find, counter- 
intuitively, that SPP scattering may increase the drift 
velocity, especially in large diameter tubes, high carrier 
density, and elevated CNT optical phonon temperatures 
due to the non-parabolicity of the bandstructure and in- 
fiuence of the higher subbands. 

Finally, it is useful to comment on the SPP scattering 
in CNTs on other possible polar gate insulators, such as 
AI2O3, AIN, Zr02, IIf02, and ZrSi04 in which the lowest 
energy phonon with an appreciable coupling is estimated 
to be hwso = SwTO\/(l + eo)/(l + eoo) = 87, 108, 37, 
24, 60 meV, respectively using the parameters in ref. [§]. 
The Frohlich couplings Eq. ^ to the low energy modes 
in these materials are similar to that in Si02. Therefore, 
we may expect that from the point of SPP scattering, 
AIN would be the best choice for the insulator, while 
Hf02 and Zr02 would give the strongest scattering at 
room temperature. 
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